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Accurate measurements of nuclear reactions of astrophysical interest within, or close to, the 
Gamow peak, show evidence of an unexpected effect attributed to the presence of atomic electrons 
in the target. The experiments need to include an effective “screening” potential to explain the 
enhancement of the cross sections at the lowest measurable energies. Despite various theoretical 
studies conducted over the past 20 years and numerous experimental measurements, a theory has not 
yet been found that can explain the cause of the exceedingly high values of the screening potential 
needed to explain the data. In this letter we show that instead of an atomic physics solution of 
the “electron screening puzzle”, the reason for the large screening potential values is in fact due to 
clusterization effects in nuclear reactions, in particular for reaction involving light nuclei. 


To understand the energy production in stars, the first 
phases of the universe and the subsequent stellar evo¬ 
lution, an accurate knowledge of nuclear reaction cross 
sections a{E) close to the Gamow energy Eg is required 
mm- Therefore, recent research in experimental nuclear 
astrophysics has triggered the development of new theo¬ 
retical methods and the introduction of new experimental 
techniques to study thermonuclear reactions at ultra-low 
energies, either directly or indirectly. In nuclear reactions 
induced by charged-particles occurring during quiescent 
burning in stars, Eg (in general of order of few keV to 
100 keV) is far below the Coulomb barrier E CB - for the 
interacting nuclei, usually of the order of few MeV. In 
particular almost all of the nuclear reactions relevant to 
solar energy generation are between charged particles and 
non-resonant reactions [2] . This implies that as energy 
is lowered the thermonuclear reactions are more depen¬ 
dent on the tunneling with an exponential decrease of 
the cross section. Therefore, their bare nucleus cross sec¬ 
tions CTh(E) drops exponentially with decreasing energy. 
For such reactions it is helpful to remove the rapid en¬ 
ergy dependence associated with the Coulomb barrier, 
by evaluating the probability of s-wave scattering off a 
point charge. The nuclear physics (including effects of 
finite nuclear size, higher partial waves, antisymmetriza- 
tion, and any atomic screening effects not otherwise ex¬ 
plicitly treated) is then isolated in the S factor, defined 
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by 


S(E) = Ecr b (E) exp [2ni](E )}, 


(1) 


where CTb(E) is the bare nucleus cross section at the cen¬ 
ter of mass energy E and exp(27r?y) is the inverse of the 
Gamow tunneling factor, which removes the dominant 
energy dependence of a b (E) due to barrier penetrability, 
and the Sommerfeld parameter ri(E) is defined as 
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It depends on the atomic numbers Zi, 7j 2 of the colliding 
nuclei and on their relative velocity v = -^/2E//i in the 
entrance channel and a = e 2 /he the fine-structure con¬ 
stant. Due to this definition, the astrophysical S factor is 
a slowly varying of E and one can extrapolate S(E ) more 
reliably from the range of energies spanned by the data to 
the lower energies characterizing the Gamow peak. (For 
more details see Refs. mm and references therein). 

The measured cross-section a s must be corrected for 
the effect of electron screening arising from the presence 
of electrons in the target atoms and, possibly, in the 
(partly)-ionized projectiles [T] [2] EHS] ■ The presence of 
electrons contribute to an enhancement of the measured 
cross-section compared to that with bare nuclei. Note 
that a similar screening effect is also present in stellar 
plasmas, where fully ionized atoms are surrounded by a 
“sea” of electrons within the so-called Debey-Huckel ra¬ 
dius, which in turn depends on conditions of plasma tem¬ 
perature and density that may vary during stellar evo¬ 
lution. Because the electron screening measured in the 
laboratory differs from the one in the plasma, it is im¬ 
portant that the measured cross sections of astrophysical 
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interest be the bare one, er;,(E), so that plasma screening 
corrections can be subsequently applied. To parameter¬ 
ize the cross section rise due to the screening effect, an 
enhancement factor fi a b{E) is usually introduced. This 
factor can be described in a simplified way by the equa¬ 
tion BEUl, 


flab(E) 
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where Ue lab ^ is the electron screening potential in labora¬ 
tory experiments. A good understanding of the electron 
screening potential C/e b is essential in order to calculate 
<Jb(E) from the experimental a s (E) using Eq. ([3]). For 
astrophysical applications it is necessary to know accu¬ 
rately the reaction rates in the stellar plasma, amounting 
to an average over the particle velocities, ( a p i(E)v ). In 
turn, the effective cross section for stellar plasma, cr p i(E), 
is connected to the bare nucleus cross section and to the 
stellar electron screening enhancement factor f p i by the 
relation 


Ppi(E) = ab{E) fpi(E). (4) 


be expressed as 

^nucleus) = Oi \i/ja) + P \i’ai’Ei) + 7 I 4>c4>d) 4- , (5) 

where (a,/3, 7 ) are spectroscopic amplitudes, | ip a) is the 
wave function of the A nucleons in a non-cluster configu¬ 
ration and I'i/’a'i/’s) represents the nuclear wave function in 
a cluster-like configuration with clusters a and B. Cluster 
configurations can alter the fusion probabilities because 
the Coulomb penetrability is reduced, as we show next. 
A simple evidence is the cross section for 6 Li + 6 Li —> 3a 
reactions at ultra-low energies which are experimentally 
found to be orders of magnitude larger than calculations 
based on barrier penetrabilities for 6 Li as non-clusterized 
spherical nuclei m and references therein). 

The basic idea of the cluster model and its relation 
to the screening puzzle is that whereas the spectroscopic 
amplitudes of cluster-like structures can be very small, 
the fusion reactions have an exponential enhancement 
for cluster-like structures. Since 6 Li can have a d + a 
cluster structure, the fusion can be enhanced because 
the Coulomb barrier for deuterons with 6 Li is suppressed. 
Due to clustering the fusion cross section can be split into 
partial cross sections in the form 


If <Jf,(E) is measured at the ultra-low energies (Gamow 
energy Eg) and f P i{E) is estimated within the frame¬ 
work of the Debye-Hiickel theory it is possible from Eq. 
© to evaluate the cr p i(E) which is the main quantity 
necessary for astrophysical applications. Unfortunately, 
direct experiments to measure the cross sections of re¬ 
actions involving light nuclides have shown that the ex¬ 
pected enhancement of the cross section at low energies 
connected to the screening effect is, in many cases, signif¬ 
icantly larger than what could be accounted for by avail¬ 
able atomic-physics models mm- This aspect deserves 
special attention ; 3] because one may have a chance to 
predict the effects of electron screening in an astrophys¬ 
ical plasma only if it is well understood under labora¬ 
tory conditions (Eqs. (l][3)). To explain the laboratory 
screening puzzle many experimental mmm and theo¬ 
retical studies have been carried out mm- In particular 
many experiments were performed to estimate the sys¬ 
tematic errors in the determination of the astrophysical 
factor. Special investigations were carried out to rule out 
errors that might be present in the extrapolation of the 
data to zero energy and in the calculations of the energy 
loss at these ultra-low energies m ei ms]. But up to 
now, theoretical studies from the point of view of atomic 
physics have not given a solution to this puzzle. This 
lack of theoretical understanding can jeopardize the sig¬ 
nificance of some values of the bare nucleus astrophysical 
factor Sb(E ) extracted from direct measurements. The 
aim of this letter is to uncover a novel approach to the so¬ 
lution of this puzzle by the introduction of nuclear struc¬ 
ture effects without questioning the well known atomic 
physics effects. The main motivation to justify the intro¬ 
duction of this new idea is explained next. 

The wave function of a nucleus in the Fock space can 


cl — CeePeu+ 6 u + C-2ePd+ 6 u + • • ■ * ( 6 ) 


where the constants Cee and C 26 include the spectro¬ 
scopic amplitudes and appropriate phase factors for 6 Li 
+ 6 Li and d + 6 Li configurations. Because 


Pe Li+6Li 0 
-Pd+ 6 Li 


(7) 


as k —> 0, the second term in Eq. ([6]) dominates at lower 
energies. Thus, even if the spectroscopic amplitudes are 
small, the cluster-like configuration enhances the cross 
section by many orders of magnitude, more than compen¬ 
sating the small configuration probability, also related to 
the preformation factor. The experimental data clearly 
shows that the fusion cross section in the 6 Li + 6 Li —> 
3a does not decrease as fast with energy as the penetra¬ 
bility for 6 Li + 6 Li channel does An additional 

effect might be responsible for an even larger enhance¬ 
ment of the cross section. When the cluster-like struc¬ 
ture in 6 Li aligns so that the two deuterons and the two 
alphas are located along a line with the two deuterons 
closer to each other right before the reaction occurs, the 
barrier for the deuterons is reduced due to its larger dis¬ 
tance to the alpha particle. Even without alignment, the 
average over all configurations is still reduced. We show 
this with a simple model based on the clusterization of 
6 Li. Dramatic effects of clusterization can be imprinted 
on the quantum tunneling probability. Consider a spher¬ 
ical coordinate system, depicted in Fig. [T] placed with 
the origin in the center of mass of the projectile with 
the z-axis along the direction of bombardment. The 6 Li 
projectile is partitioned into two clusters, a deuteron (1 
in the following) and an alpha particle (indicated with 
2), with centers of mass at a distance d from each other 
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FIG. 1. Top figure: Coordinate system used in calcula¬ 
tions. Plot - Upper panel: Ion-ion potentials for the diclus¬ 
ter systems with three different orientations as a function of 
c.m. distance to show the change in barrier height and posi¬ 
tion. Plot - Lower panel: Penetration probability of averaged 
dicluster-discluster system (solid) compared with sphere-on- 
sphere (dashed). 


on the 0 -axis, in such a way that the deuteron is on the 
positive side and therefore at distances +2d/3 and —d /3 
respectively. The target center of mass is found at an an¬ 
gle 9 with respect to the projectile’s reference frame. The 
azimuthal angle of the projectile does not play any role, 
because it amounts to a rotation around the 0 -axis. The 
orientation of the generic inter-cluster axis of the target 
(here an identical 6 Li system) has angles 9' and </>'. 

The position, height and width of the barrier depend 
on the details of the potentials, sum of Coulomb and 
nuclear potentials, between each pair of clusters and it 


TABLE I. The experimental values of the electron screening 
potentials, Ut xv , and theoretical adiabatic limits, lJ e adlzrn _ 



Reaction 

jj-adlim 

U e e xv 

Note 

Ref. 



(eV) 

(eV) 



[1] 

2 H(d, ffH 

14 

19.1±3.4 


m\m 

[21 

3 He(d,p) 4 He 

65 

109±9 

D 2 gas target 

m 

131 

3 He(d,p) 4 He 

120 

219±7 


m 

[4] 

3 He( 3 He,2p) 4 Hc 240 

305±90 

compilation 

m 

[51 

°Li(d,a) 4 He 

175 

330T120 

H gas target 

m 

161 

6 Li(d,a) 4 He 

175 

330±49 


Bung 

[7] 

6 Li(p,a) 3 Hc 

175 

440T150 

H gas target 

m 

181 

6 Li(p,a) 3 Hc 

175 

355±67 


m Isa nn 

[9] 

7 Li(p,a) 4 He 

175 

300T160 

H gas target 

m 

[10] 

7 Li(p,a) 4 Hc 

175 

363±52 


on nn Ha 

[HI 

a Be(p,a 0 ) b Li 

240 

788±70 


HIES] 

[12] 

10 B(p,« 0 ) 7 

340 

376±75 


[MEZ] 

[13] 

11 B(p,a 0 ) 8 Be 

340 

447±67 


EH EH] 


depends on vector distances r,y as 


Vu, t(r, 9 , 9\ 0') = E ( + Mr«) 

ij=i V / 


( 8 ) 


where the function parametrically depends also on the 
relative inter-cluster distance d , that we keep constant 
(and equal to ~ 3.85 fm obtained from the cluster model 
formula (A.4) of Ref. [13] and data on radii [H]). The 
next level of refining of the model would include the 
weighting with the internal relative motion wavefunction 
of 6 Li. The upper panel of Fig. [l] shows the significant 
changes on the barrier height and position induced by 
relative rotations of the di-cluster orientation axes, while 
keeping 9' = ft = 0. 

For energies below the barrier, the tunneling proba¬ 
bility can be calculated in the WKB approximation as 
P = e ~ 2G , where the Gamow tunneling factor is given 
by 


G{E,9,6’A ') 



(9) 

The angle-averaged penetration probability as a func¬ 
tion of bombarding energy is displayed in the lower panel 
of Fig. |T] (solid line) and compared with the analogous 
calculations for sphere-on-sphere (optical potential from 
Ref. [15]). 

It is clear from this simple analysis that the probabil¬ 
ity is very much enhanced in the dicluster-dicluster fusion 
model with respect to the sphere-on-sphere model. Thus, 
if the spectroscopic amplitudes in Eq. ([5]) are known, 
the coefficients in the partial cross sections of Eq. ([6]) 
will also be known and the total fusion cross section will 
certainly display the enhancement effects due to cluster¬ 
ization. The problem of calculating those amplitudes is 
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FIG. 2. Ratio of the experimental electron screening poten¬ 
tial Ut xp and the theoretical adiabatic limit of the electron 
screening potential as function of the main reaction 

present in the literature. The vertical bars are the total un¬ 
certainties of the measurements reported in literature. The 
numbers in brackets correspond to those in Table I. 


a very difficult one, not within the scope of this article. 
It requires a theory beyond the naive shell model, which 
treats nuclei as a collection of nucleons. On the other 
hand, cluster models rely on the knowledge of preforma¬ 
tion factors. In this respect, ab-initio models are quite 
promising (see, e.g., Ref. [29]) but the inclusion of cor¬ 
relations including clusters in bound states has shown 
to be quite challenging and one does not seem to have 
reached the stage of properly assessing the values of the 
spectroscopic amplitudes for each cluster configuration. 
The clustering effect we propose as a candidate to explain 
the electron screening puzzle is somewhat related to the 
Oppenheimer-Phillips effect nanu, which is due to the 
polarization of the deuteron in the Coulomb field of the 
target nucleus in deuteron induced reactions. We have 
shown that even without polarization the fusion of light 
cluster-like nuclei can acquire enhanced tunneling when 
averaged over all geometric configurations. 

In Table I and figure 2 we show typical cases of reac¬ 
tions at ultra-low energies where clusterization fusion en¬ 
hancements might be have been observed: the first is for 
the case of Z = 1 nuclei reacting with nuclei which do not 
present an evident nuclear cluster structure, the second 
is for the case of cluster-like nuclei. The main conclusion 
drawn from Table 1 is that there is a clear correlation 
between the cluster structure of nuclei involved in reac¬ 
tions at ultra-low energies and the discrepancy between 
the value of the upper limit (adiabatic approximation) 
of the screening potential, and its experimental 

value, U e exp \ The disagreement increases as the clus¬ 
ter structure is more pronounced (larger cluster spectro¬ 
scopic factor). In particular, Table 1 displays the fol¬ 
lowing evidences in favor of a nuclear structure solution 
for the “electron screening” puzzle for the thermonuclear 
reactions of astrophysical interest: 


(i) - In all the cases of reactions with cluster-like nuclei 

with one electron (Z = 1 ) the experimental elec¬ 
tron screening potentials are in agreement, within 
the experimental errors, with the upper theoretical 
limit due to atomic energy balance [ 1 ] (examples 
are the cases of d + d and d + p reactions). 

(ii) - For reactions with cluster-like nuclei with small 

electron number (Z = 3) and mass numbers 6,7 
(examples are p + 6,1 Li) if we consider the central 
values of the experimental screening potential we 
observe that these values are more than a factor 
1.5 times higher with respect the \J e adhm . 

(iii) - Reactions with cluster nuclei with electron number 

Z = 4 — 5 also show a disagreement between the 
experimental and theoretical upper limit based on 
the energy balance in the adiabatic approximation. 
The discrepancy increases as the cluster structure 
of the interacting nuclei is more evident (examples 
are the reactions p + 9 Be and p + 10,11 B). 

If the solution of the “electronic screening” puzzle would 
be related to the effects of atomic nature the item (ii) 
of the list above should have the same value of exper¬ 
imental screening potential for all three reactions (p + 
6 Li, d, p + 6 Li and 7 Li reactions) because of the isotopic 
invariance. In fact, in the case of the reaction p + 7 Li 
there is a deviation of the central value of experimental 
screening potential by about 250 eV, while for other cases 
there is a deviation of 165 eV (d + 6 Li) and 180 eV (p + 
6 Li). If we consider the case of the reactions in the group 
(iii) we find that the differences found with the atomic 
screening predictions are 90 eV for p + 10 B, and 90 — 130 
eV for the reaction n B + p. For these latter cases, in 
order to draw more definitive conclusions, new measure¬ 
ments are needed with larger precision to reduce the mar¬ 
gin of uncertainty. In any case, these results are not in 
agreement with the atomic description of the screening of 
thermonuclear reactions in the laboratory. The main idea 
introduced in this work is that at the very low energies of 
nuclear reactions in stellar environments, the condition 
set by Eq. 0 may occur due to the presence of clusters in 
the interacting nuclei. The Coulomb penetrability is sup¬ 
pressed with decreasing energy of interacting nuclei due 
to clustering and polarization. Therefore, the absolute 
values of astrophysical factors obtained through extrap¬ 
olation should be reviewed to include not only atomic 
physics but also nuclear physics effects to correct for the 
increase of the astrophysical S-factor as the energy de¬ 
creases. Indeed, in Eq. (|T|) it is implicitly assumed that 
the wave function for the relative motion of the nuclei 
is expressed only by the first term of Eq. 0 (a = !)• 
Therefore the Gamow factor exp(27r?7) is calculated tak¬ 
ing into account that in Eq. 0 no cluster structures 
exist, i.e. /3 = 7 = 0. Only in the situation that clusters 
can be formed and some sort of polarization occurs, the 
Gamow factor can compensate for the drastic suppression 
of the cross section with decreasing energy in many of the 
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reactions of astrophysical relevance. These conclusions, 
to be confirmed by further more precise experiments, will 
lead to a critical review of the actual values of the elec¬ 
tronic screening potentials. This problem also appears to 
exist with direct experiments at higher energies. From 
the considerations expressed in the present work we can 
state that the discrepancy between the experimental elec¬ 
tron screening potential values and the upper theoretical 
values (adiabatic limit) may be linked to nuclear struc¬ 
ture effects and not to hitherto unknown and speculative 
processes in atomic physics. We propose that new the¬ 
oretical and experimental studies in the field of nuclear 
astrophysics at very low energies should be carried out. 
In particular, a more comprehensive theoretical reaction 


method that takes into account polarization and align¬ 
ment of cluster-like nuclei should be pursued. New and 
more precise measurements to confirm this theory should 
be carried out concomitantly. The nuclear reactions in¬ 
volving 6 Li and 7 Li, such as the 6 Li + 6 Li, 'Li + 7 Li, 
7 Li + 6 Li , 9 Be + 3 He, 9 Be + 'Li should be of particular 
interest to prove the relevance of such additional nuclear 
structure effects in thermonuclear reactions. 
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